Search for quantum dimer phases and transitions in a frustrated spin ladder 
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A two- leg spin- 1/2 ladder with diagonal interactions is investigated numerically. We focus our 
attention on the possibility of columnar dimer phase, which was recently predicted based on a 
reformulated bosonization theory. By using density matrix renormalization group technique and 
exact diagonalization method, we calculate columnar dimer order parameter, spin correlation on a 
rung, string order parameters, and scaled excitation gaps. Carefully using various finite-size scaling 
techniques, our results show no support for the existence of columnar dimer phase in the spin ladder 
under consideration. 

PACS numbers: 75.10.Pq, 75.10.Jm, 75.30.Kz, 75.40.Cx 



Much effort has been devoted to understanding the ef- 
fects of competing interactions on quasi-one-dimensional 
systems. The possibility of unconventionally ordered 
phases has been the focus of interest. For example, one- 
dimensional extended Hubbard model (EHM) with near- 
est neighbour repulsion V in addition to on-site repulsion 
U has been investigated intensivelyi 1 i 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 i 10 T 11 i 12 
It had been considered for a long time that the ground- 
state phase diagram at half filling has only two phases, 
the spin-density-wave (SDW) and the charge-density- 
wave (CDW) states! 1 ' 2 ' 3 ! 4 ! 5 ! 6 Moreover, the order of the 
phase transition at U ~ 2V can change from continu- 
ous to first order at a tricritical point, which was specu- 
lated to exist in the intermediate coupling regimei^* 3 ^ 
Later, by using the level-spectroscopy method, Naka- 
mura pointed out that there exists also a novel sponta- 
neously dimerized phase, the so-called bond-order-wave 
(BOW) phase, in a narrow strip between the SDW and 
the CDW phases in the weak coupling regionJ While 
this phase is absent in standard one-loop g-ology and 
bosonization calculations) 1 ^ its existence was supported 
by a reformulated one, 8 where higher-order terms were 
included. The appearance of the BOW phase was sub- 
sequently confirmed by quantum Monte Carlo simula- 
tionsj2ii& density matrix renormalization group (DMRG) 
method^* and functional renormalization group calcula- 
tions^ 

Interestingly, same story may happen also in an anti- 
ferromagnetic two-leg spin ladder with diagonal frustra- 
tions. The Hamiltonian of this model reads as follows, 



H — J S a ,i • S Q i + i + J± Si,, ■ S2,i 

a— 1,2 i i 

(Si,j • Sa,i+i + Sj.,4+1 • Sa,i) , (1) 



where S a j denotes a spin-1/2 operator at site i of 
the a-th leg. J and J± are the exchange couplings 
on legs and rungs, respectively. We set J = 1 
hereafter. J x denotes the next-nearest-neighbor inter- 
chain coupling. This model has been investigated for 
decadesiiiii 4 *i£*i 6 ^i£ii2ii2i2S and considerable amount of 



knowledge has been accumulated. It has been believed 
that the ground state phase diagram consists of only 
two phases, the rung-singlet (RS) phase and the so- 
called Haldane phase. The earlier bosonization study 
predicted a direct transition between these two phases 
at J± = 2J X . 14 Various numerical calculations, such 
as series expansions, 18 DMRGfiS and level-spectroscopy 
method, 20 showed the phase boundary to be shifted away 
from the J± = 2 J x line, with a larger RS phase. More- 
over, the phase transition seems to be of second order 
at weak interchain couplings and becomes first order at 
stronger couplings. 18 - 19 Recently, it was suggested that 
there exists an intermediate, spontaneously dimerized 
phase, the so-called columnar dimer (DC) phase, lying 
between the Haldane and the RS phases at weak inter- 
chain couplings. 21 This remarkable proposal is based on 
a reformulated weak-coupling field theory, which is simi- 
lar in spirit to that in Ref . with success in the study of 
one-dimensional EHM. According to Ref. 0, for a given 
J x , the DC phase occurs within a narrow but finite re- 
gion (J±)c,t < J± < (J±)c,s in the phase diagram. Here 
(J±)c,t = 2J X -5J 2 /7r 2 and (J±) c ,s = 2J X - J 2 -/ 71 " 2 are 
two distinct critical points, given by vanishing mass gaps 
in the spin-triplet and the spin-singlet sectors, repsec- 
tively. As mentioned before, this intermediate DC phase 
was not found in previous numerical calculationsii&iSiSSi 
Thus the existence of this spontaneously dimerized DC 
phase is surprising and calls for thorough theoretical 
studies. 

In this paper, we try to find numerical evidence for 
the possibility of the DC phase in a two-leg spin ladder 
of Eq. JIJ in the vicinity of J± — 2 J x . Here we take J x = 
0.2m£ It is smaller than the value J x = 0.287, where the 
transition changes to be first order J£ Thus the proposed 
region for the DC phase becomes 0.38 < J± < 0.396. 
The values of J x and J± under consideration should 
be small enough in accordance with the weak-coupling 
field theory in Ref. |U Both the DMRG technique 2 ^ 
under open boundary conditions up to L = 400 rungs 
and the Lanczos exact diagonalization method with pe- 
riodic boundary conditions up to L = 16 are used. In our 
DMRG calculations, 500 states per block are kept, and 
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FIG. 1: (Color online) (a) Size dependence of the DC order 
parameter D for various J± with J x = 0.2. (b) In D versus L 
for various values of J± with J x = 0.2. Labels for various J± 
are the same as those in (a). Inset shows the inverse of the 
correlation length £ for various values of J±. 



the truncation error is of the order of 10~ 7 . To demon- 
strate the possibility of the DC phase, the most direct 
way is to show the corresponding order parameter be- 
ing nonzero in the proposed region of the phase diagram. 
As far as we know, this order parameter has not yet been 
measured for the present model. We note that the critical 
point (J±) c ,t in the spin-triplet sector is consistent with 
the phase boundary found previously P i&iiL2 n i Therefore, 
another support of the proposal in Ref. 21 is to show 
the existence of another critical point at (J±) c .s with 
vanishing mass gap in the spin-singlet sector. Carefully 
using various finite-size scaling skills for the DC order 
parameter and other physical quantities defined below, 
our results fail to show the existence of the DC phase, 
but instead they indicate a direct transition between the 
RS and the Haldane phases. 

By using the DMRG technique under open boundary 
conditions, we first analyze the DC order parameter. As 
shown in Ref. l25l due to the presence of open ends, weak 
dimerization profiles can be induced near the boundaries. 
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FIG. 2: (Color online) Spin correlation on the (iV/2)-th rung 
as a function of J± for various sizes L with J x = 0.2. Inset 
shows the scaling behavior of the crossing points. 



In order to reduce the boundary effect, the DC order 
parameter is calculated by the difference of local spin 
correlation on legs, 




where only the bonds in the middle of finite open lad- 
der with length L are considered (i.e., i = L/2). (• • •) 
means the ground-state expectation value. The DC or- 
der parameter in the thermodynamic limit is then = 
limL—nx D. Our results of the DC order parameter D 
for various J± with J x = 0.2 are shown in Fig. ^ ( a )- 
We find that D always decreases to zero in the thermo- 
dynamic limit even for J± — 0.39, which lies within the 
suggested region for the DC phase. Fig. ^ (b) shows 
lnl? versus L for various values of J±. It is found that, 
for L being large, D « cexp(— L/£), where c is a con- 
stant and £ is a kind of correlation lengths. Moreover, 
as shown in the inset of Fig. ^ (b), l/£ oc \ J± — («7j_) c | 
with (J_l) c ~ 0.38, which agrees with the value of the 
proposed critical point (Jj_) c ,t = 0.38 in the spin-triplet 
sector. This indicates that the long-range DC phase may 
appear only at this phase boundary, rather than for a 
finite region in the phase diagram. Besides, it shows no 
evidence for the additional second-order quantum phase 
transition at J± = (J±) c ,s = 0.396, since the correlation 
length £ diverges only at a single critical point (J±_) c .t, 
rather than at two points. 

To be sure if we miss the critical point (Jj_) c ,s in the 
above analysis, a finite-size crossing method 26 is used, 
which is applicable to detect the critical points of second- 
order quantum phase transitions. It is noted that one can 
always decompose the Hamiltonian into two parts, i.e., 
H = Tio+gV, and consider the transition being driven by 
the parameter g. Based on the finite-size scaling analysis, 
it is shown that the curves of the mean value O = (V) / L 
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FIG. 3: (Color online) String order parameters for ladders 
with various sizes L as functions of J± with J x = 0.2. Inset 
shows the L _1 scaling behavior of the crossing points. 



at two successive values of size L as a function of g will 
cross at a single point g* L near each critical point g c S^ 
The value of the critical point g c can be found numer- 
ically by extrapolating the sequence g* L to L — » oo. In 
the present case, we take the driving parameter as Jj_. 
Thus the corresponding transition-driving term becomes 
V = J2i S M ' s 2,i, which gives O = (Si,, ■ S 2 ,i) by trans- 
lational invariance if periodic boundary conditions are 
used. In case of finite open ladders, to avoid bound- 
ary effects, sites in the middle of ladders (i.e., i = L/2) 
are used. In Fig. [21 we plot the curves of (Si,j • S2,i) 
versus Jj_ for various sizes L, which are calculated by 
the DMRG technique. It is found that there is only one 
crossing point Jj_(L) at L = (L\ + L,2)/2 for the curves 
at two subsequent sizes L\ and L^. This indicates that 
there is only one, but not two, phase transition. Our con- 
clusion is consistent with the picture provided by previ- 
ous investigationSfiiiiiiSiiSiiLi2ii2i2& where the complete 
phase diagram consists of only two phases with a single 
phase boundary. The scaling behavior of the crossing 
points is shown in the inset of Fig. [21 It is found that 
the crossing points converge to the value (J±) c — 0.378 
for the critical point. Our finding agrees well with that 
obtained by previous DMRG calculation— and is consis- 
tent with the predicted value for (J±) c ,t- 21 Again, our 
results provide no support for the existence of another 
critical point at (Jj_) c ,s in the spin-singlet sector. 

While our results show that there is no DC phase and 
there are only two phases in the complete phase diagram, 
the nature of these two phases has not yet been explored 
in the present study. According to previous investiga- 
tions, these two phases should be the Haldane and the 
RS phases, and they can be identified by two distinct 
string order parameters O dd and O cven ii^iSiiL21 These 



two string order parameters are given by 

Op = - ^ ^%exp ^tt £ S Z P ^ S Z P ^ , (3) 

where P — odd, even. The composite spin operators are 
defined as Sl dA i = Sf^ + Sf ,» an d S* vcn ^ — + of,i+i- 
Because of spin isotropy, we calculate the string order pa- 
rameters for the z-component spins only. In case of finite 
ladders, it turns out that the intersection of the curves of 
two distinct string order parameters implies the critical 
point ii£ In order to reduce the undesirable boundary ef- 
fects, we fix site j in Eq. at the center of the chain and 
let site i = 20 in our calculations. Our DMRG results 
of Oodd and Ooven for various sizes L as functions of J± 
with J x = 0.2 are shown in Fig. [31 For smaller J±, one 
has O dd 7^ and O e vcn = 0, which implies the Haldane 
phase; while dd = and O CV cn ^ for larger J± , which 
implies the RS phase. The finite size scaling procedure is 
used to determine the thermodynamic limit of the value 
of the critical point. As shown in the inset of Fig. [21 the 
crossing points converge to the value of {J±) c — 0.378, 
which agrees quite well with that found in the inset of 
Fig. El 

Finally, we provide a further examination of the pos- 
sibility of another critical point (J±) c ,s with vanishing 
gap in the spin-singlet sector. According to the phe- 
nomenological renormalization-group (PRG) method^ 
second-order phase transition points can be determined 
by the crossing points of the curves of the scaled gaps 
LAE V at two successive sizes L and L + 2, where AE V 
denote the excitation gaps in the spin-singlet and the 
spin-triplet sectors for v = s and t, respectively. Here 
the gaps are calculated by using exact diagonalization 
method with periodic boundary conditions up to L = 16. 
In the present case, the ground state is unique for any 
J x and J±, and it has total spin S = and momen- 
tum k = 0. The spin-triplet (-singlet) excitation gap 
AEt {AE S ) is determined by the energy difference be- 
tween the ground state and the lowest level with total 
spin S = 1 (with total spin 5 = and k — n). Our 
results for the scaled gaps LAE V as functions of J± for 
various sizes L with J x = 0.2 are exhibited in Fig. 01 For 
both cases of LAE S and LAE t , there are two crossing 
points of the curves at subsequent sizes. However, it im- 
plies only one, but not two, critical point, because the ex- 
trapolation of the left and the right crossing points tend 
to a single value in the thermodynamical limit as shown 
in the insets of Fig. 01 The limiting values for the spin- 
singlet and the spin-triplet sectors are almost the same, 
and both give (J_l) c — 0.38. Our findings indicate that 
only one phase transition at (Jj_) c — 0.38 occurs in the 
present system, and then the spin-singlet and the spin- 
triplet gaps vanish simultaneously. We note that simple 
linear extrapolations from data of systems of small sizes 
may be somewhat dangerous. Nevertheless, because our 
results based on PRG are consistent with that obtained 
by the above DMRG analysis and that given by earlier 
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DMRG calculation^ they may lead to true physics in 
the thermodynamical limit. 



In summary, we study numerically a two-leg spin 
ladder with diagonal frustrations for weak interchain 
couplings. All our results indicate that the transition 
from the RS to the Haldane phases is direct, with- 
out any phase in between. This conclusion is consis- 
tent with the picture obtained by all previous investiga- 



tions, 



13,14,15,16.17,18.19.20 



except that proposed in Ref . |21| 



based on a reformulated weak-coupling field theory. It is 
not clear why a reformulated bosonization analysis works 
for one-dimensional EHM, but could fail for the present 
two-leg spin ladder. Further theoretical investigations 
are necessary to clarify this issue. 



FIG. 4: (Color online) Scaled gaps (a) LAE 3 for spin-singlet 
excitation and (b) LAE t for spin-triplet excitation as func- 
tions of J± for various sizes L with J x =0.2. Insets show the 
L~ 2 scaling behavior of the left and the right crossing points 
represented by open and solid circles, respectively. Dotted 
lines are guides to the eye. They are straight lines fitted to 
the last few points. 
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